Spin-charge separation in strongly interacting finite ladder rings 
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We study the conductance through Aharonov-Bohm finite ladder rings with strongly interacting 
electrons, modelled by the prototypical t-J model. For a wide range of parameters we observe 
characteristic dips in the conductance as a function of magnetic flux, predicted so far only in chains 
which are a signature of spin and charge separation. These results open the possibility of observing 
this peculiar many-body phenomenon in anisotropic ladder systems and in real nanoscopic devices. 
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I. INTRODUCTION 

The phenomenon of the fractionalization of an elec- 
tron into its spin and charge degrees of freedom was pre- 
dicted tlieoretically for one dimensional (ID) strongly 
interacting systems in the framework of the Luttinger 
Liquid (LL) theoryJ^i^ Continuing progress in fabrica- 
tion techniques, and the discovery of new materials of 
quasi- ID electronic character, have led in the last decade 
to a variety of experiments which seek evidence of spin 
charge separation (SCS) such as the observation of non- 
universal power-law I-V characteristics?^ the search for 
characteristic dispersive features by angle-resolved pho- 
toemission spectroscopy (ARPES),^ the violation of the 
Wiedemann-Franz law, — and the analysis of spin and 
charge conductivities i^ii Among the candidate materi- 
als to present SCS^ we can mention the organic Bech- 
gaard and Fabre salts, molybdenum bronzes and chalco- 
genides,— cuprate chain and ladder compounds^, later- 
ally confined two-dimensional electron gases, cleaved- 
edge overgrowth systems^" and also carbon nanotube sys- 
tems J^ii^ 

From the theoretical point of view, several ways for 
detecting and visualizing SCS were proposed. Direct cal- 
culations of the real-time evolution of electronic wave 
packets in Hubbard rings revealed that the spin and 
charge densities dispersed with different velocities as an 
immediate consequence of SCS.^^'^** The analysis of the 
electronic transmission through Aharonov-Bohm (AB) 
rings -'^^i-'^^i-'^'^ described by a LL presented striking fea- 
tures characteristic of SCS, where the flux-dependence of 
the transmission was found to show new structures ap- 
pearing at fractional flux values in addition to the non- 
interacting flux quantum periodicity $0 = hc/e. Ini^ 
these fractions were determined by the ratio between the 
spin and charge velocities Vg/vc- In their interpretation 
the dips arise because transmission requires the separated 
spin and charge degrees of freedom of an injected electron 
to recombine at the drain lead after traveling through the 
ring a different number of turns in the presence of the AB 
flux. However, recent results which go beyond the sin- 
gle pole approximation used there, suggest that this idea 
is too simple and claim that the number of dips is not 
determined approximately by v^/vs but by vj/v^, where 



vj is the current velocity.!^ The results, however, agree 
for small integer values of p and q, for Vg/vc — p/q- Re- 
cent numerical calculations of the transmittance through 
finite AB rings described by the t-J model show clear 
dips at the fluxes that correspond to the ratio Vsjvc^ 
As we explain below, the discrepancy arises due to the 
finitencss of the system. 

In spite of the clear indications of the existence of spin- 
charge separation in ID interacting systems and its ab- 
sence in three dimensions where the Fermi Liquid theory 
is valid, there is no final word for two dimensions. The 
non- Fermi-liquid normal state properties of high temper- 
ature superconductors have led to attempts to trace their 
origin in the possible realization of SCS in strongly cor- 
related electron systems in 2Dii^ 

In this paper we explore the possibility of the existence 
of SCS in ladders, as a first step towards two dimen- 
sions. We analyze the conductance through rings formed 
by two-leg ladder systems described by the t-J model as 
a prototype of interacting systems. For certain parame- 
ters we find, indeed, clear dips at fractional values of the 
magnetic flux which we can interpret as fingerprints of 
charge and spin separation due to the difference in the 
charge and spin velocities. 

II. THE MODEL 

Our model Hamiltonian reads H = i?ieads + ^^link + 
Hring, (Fig- 1), where i/icads describes free electrons in 
the left and right leads, 

i^link = -t'^ia-l,a(^Oi,<T + aI,o-CL/2i,<T + H.c.) (1) 
(J 

describes the exchange of quasiparticles between the 
leads (fli^cr) and particular sites of leg 1 (ci^^o-), and 

Hnng = -eVg J2 - i||(4,,c,,+i,.e-^*/^ -f H.c.) 

- i-LE(cL-C^.,a + H.c.) + i/i„t (2) 

describes the interacting electron system. The fermionic 
operators c]^ ^ create an electron at site i — 1,L of leg 



FIG. 1: Schematic representation of the system. 



Z = 1, 2 with spin a. The AB ring has L rungs, is threaded 
by a flux (j) {(p — 27r"I>/$o)i and is subjected to an apphed 
gate voltage Vg. 

As we will consider the t-J model for the ladder, the 
interacting part of the Hamiltonian reads: 



Hi, 



(3) 



1,1 



where S^, = cJ,Q'''a/3Cii/3 is the spin at site i and leg 
I and no double occupancy is allowed. 

When the ground state is non-degenerate, the zero 
temperature transmission from left to right can be cal- 
culated to second order in t' by means of the retarded 
Green function for the isolated ladder ring between sites 
i and j: G^j{u!), for an incident particle with energy uj 
and momentum zfcfci ^^i^^ 



T{uj,Vg,> 
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where e{uj) — t'^G^_^ Oii^)' the effective hopping across 
the ring is i{Lu) = t'^G^^ j^^^^ (w), and uj = —2t cos k is the 
tight-binding dispersion relation for free electrons in the 
leads. This equation is in fact exact for a non-interacting 
system; with interactions on the ring it serves as an ap- 
proximation in the tunneling limit t' /t <ti 1^^^^^ for a 
non-degenerate ground state. For an odd number of elec- 
trons, the ground state is Kramers degenerate for a sys- 
tem with time reversal symmetry and the equation ceases 
to be validii^iSS So we assume that the ensuing Kondo 
effect is destroyed by temperature or magnetic HeldM^ 
The conductance is G = {2e^/h)T. 

We calculate T(LO,Vg,(j)) by numerically diagonalizing 
the isolated interacting ring in the presence of a magnetic 
field with L rungs and N electrons in the ground state 
to obtain the Green functions which appear in Eq.(|4]), 
fixing the chemical potential of the non-interacting leads 
to zero {uj = 0). By varying Vg, T{0,Vg,(j)) presents 
narrow peaks with a width proportional to {t')'^ at gate 
voltages which correspond to the excitation energies of 
the system with — 1 particles.!^'-- The transmittance 
is obtained by integrating the spectra over a small energy 
window at the Fermi energyJ^ii^iii 

III. RESULTS 

In order to study the robustness of the spin-charge sep- 
aration in the presence of a second chain, we first show 




FIG. 2: Weakly coupled chains: Transmittance as a function 
of flux for the anisotropic t-J model with J — 0, — 0.1, 
L = 6 rungs, t' = 0.05t|| , several values of t± and TV = 6 (top) 
and N — 8 (bottom) electrons in the ground state. 



results for weakly coupled chains {t± ^ ty) and Ji = 
(i II), for which we know that in one chain there 
is complete SCSj^^'^'^^ In Fig. 2 we show the results 
for several small values of t± and in fact, observe clear 
dips at certain fractional values of the magnetic flux (the 
abrupt jumps correspond to other level crossings). This 
is the first evidence, up to our knowledge of charge-spin 
separation in finite systems with more than one chain. 

To understand the position of the dips we resort to 
the expression obtained in ID for J = {i =_L, 
Considering a non-degenerate ground state containing 
A^ = A^e + 1 particles and analyzing the part of the Green 
function that enters the transmittance when a particle 
is destroyed, it is shown that the dips occur when two 
intermediate states cross at a given flux and interfere de- 
structively. These particular fluxes depend on the spin 
quantum numbers, and are located at 



(j)d = 7r(2n+ 1)/A^e 



(5) 



with n integer. If the integration energy window includes 
these levels, a dip in the conductance arises. 

For the ladder with t± = and a total even number 
of electrons A" in the ground state, the lowest-lying state 
has A^/2 electrons in each leg. As we are calculating the 
transmittance through one leg only, and the intermediate 
state has one particle less, from the condition for (pd with 

In 



A'e = N/2 — \, one expects to see dips ai (pd = 



- ~ 2n+l 
N/2-1 ■ 

Fig. 2 we see that this is the case since for the top figure 
there will be A'e 4- 1 = 3 electrons in each leg, leading 
to a dip at (p — tt/2 and for the bottom figure there will 
be 4 electrons in each leg leading to a dip at (/> = 7r/3. 
When the "second" dimension is turned on and ij^ 7^ 0, 
we find that the dips remain and are quite robust, even 
for values of t±/t\\ as high as 0.1. 

These are so far the results for a weak interchain cou- 
pling. Another atracting case is the one with a large 
coupling between the legs, ij_ 3> . In this limit and for 
the noninteracting case, the bands corresponding to the 
bonding and antibonding states of each rung are very far 
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FIG. 3: Strongly coupled chains: Flux-dependent transmit- 
tance for L = 6 rungs, J = 0, t\\ — 0.1. Top: t' — 0.05f|| and 
N = 4l electrons. Bottom: t' = 0.09tii and = 6 electrons. 



apart and one might expect the reappearance of SCS. In 
fact, this is the case, as can be seen in Fig. 3, where we 
plot the transmittance for a ladder with several values of 
tj^ and fillings. Now the total number of electrons in the 
lower band corresponds to the total filling N (for a less 
than half filled band) and the transmittance will involve 
Ne = N —\ electrons. Hence, if SCS exists, the dips will 
be found at fluxes <f)d — ""T^^f- In this figure we find 
that for large values of t^/t\\ , the dips correspond indeed 
to these fluxes. For smaller values of t±, we find a shift 
in the location of the minima and sometimes a splitting 
of the dips. 

It is interesting to visualize the behaviour of the dips 
when changing the parameters from weakly interacting 
chains to the strong coupHng case (small to big t_\_/t\\). 
In Fig[4] we collect the data for = 6 particles, where 
we see that for ^ 1 only one dip shows up (re- 

flecting the bahaviour of half the number of particles in 
each band corresponding mainly to each leg of the lad- 
der as in Fig. 2). On the contrary, for t_\_/t\\ ^ 1 all 
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FIG. 4: Transition from weakly-coupled to strongly-coupled 
chains: Flux-dependent transmittance for several relations of 
inter to intrachain hoppings tx/ty for L = 6 rungs and N = 6 
particles. 
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FIG. 5: Transmittance through a bonding channel and com- 
parison with Hcff for the same system as in Fig. 3. 



N particles belong to the relevant lower (bonding) band 
and two main dips arise at the positions corresponding 
to this strongly coupled case as in Fig. 3. 



IV. EFFECTIVE MODEL FOR STRONGLY 
COUPLED CHAINS 

In order to study the strongly coupled case and un- 
derstand the transition towards the limit of isotropic ex- 
change, we have mapped our ladder Hamiltonian i?ring to 
an effective model in the subspace of the bonding states 
of each rung, using degenerate perturbation theory up to 
second order in ^ The model is valid for energies lower 
than t± and a less than half-filled system: 

Has = -t||^(cLQ+i,. + H.c.)-f J^(S,-S,+i-l/4) 
+t" ^(cI+2^,Q,(SrS,+i - 1/4) + H.c), (6) 

where J = J\\/2 + 2tj/{tj_-3Jj_/A), t" = tjj/(<_L -3Jj^/4) 
and Cia — -^{ci^.cr + Ci^,a), the bonding operator. 

The second dimension of the original ladder is reflected 
by the second nearest-neighbour term which tends to de- 
stroy the dips. In Fig. 5 we show the transmittance of a 
particle through a bonding channel and compare it with 
that of iJeff, finding excellent agreement for the strongly 
coupled case. As t±/t\\ diminishes, the curves start to 
differ as Hcs loses its validity. Comparing with Fig. 3, 
where a particle is injected to one site only instead of 
into a bonding state, we find that the effective model is 
also valid for this case, as long as t± ^ iy, since the 
antibonding band is shifted to upper energies. 

So far we have presented results for Jy = J± = 0. Fi- 
nite interactions introduce an extra spin shuffling in the 
system, mixing the spin wave quantum numbers. In spite 
of the fact that the conditions that lead to Eq. (5) (based 
on the J = limit) are no longer valid, dips are still ob- 
served for small values of the interactions in the strongly 
coupled case. However, the effect of the interaction is to 
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FIG. 6: Transmittance through a bonding channel (lines) 
and comparison with Hcs (symbols) for finite interactions Jy 
and Jj_ for t_i_ — 50t|| , L — 6 rungs and TV = 6 electrons. 

reduce the depth of the dips and shift their position. In 
FiglB] we show results for the 6-rung ladder with = 6 
particles in the ground state and several values of the 
interactions. Taking into account the fact that the J's 
are obtained perturbatively from the large-?/ Hubbard 
model, we keep their relation as J±/J\\ = f^/^- Two 
observations can be made: i) The dips are still present 
for finite J's. ii) However, as for the J = case where the 
dips were affected by the interchain hopping parameter, 
in this case we also find shifts and reductions in the their 
depth caused by the interactions. A similar behaviour 



occurs for weakly interacting chains. We also find that 
the effective model fits quite well the results in the ladder 
for finite J± and J|| , and its range of validity extends to 
appreciable values of the interaction parameters. 



V. CONCLUSIONS 

In summary, we have found that the dips in the conduc- 
tance, predicted to appear in strongly correlated chains 
as a consequence of spin-charge separation, are robust in 
the presence of a second transmission channel modelled 
by a ladder system in the anisotropic limit. For a wide 
range of parameters, in particular for weak and strong 
hopping couplings across the rungs t±, the dips remain. 
However, their position differ from the predictions stem- 
ming from the exactly solvable case of the Hubbard chain 
with infinite U (or t-J model with J = 0).^^ For inter- 
mediate values of t± the dips disappear. The signatures 
of spin-charge separation are also robust for finite, al- 
beit small, values of spin-spin interactions. This is the 
first time in which signatures of spin-charge separation 
are observed in interacting finite systems with more than 
one chain, thus opening the possibility of measuring this 
peculiar phenomenon in real nanoscopic systems. 
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